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ABSTRACT 

We measure the projected spatial correlation function w p (r p ) from a large sample combining GALEX 
ultraviolet imaging with the SDSS spectroscopic sample. We study the dependence of the clustering 
strength for samples selected on (NUV — r) a f, s color, specific star formation rate (SSFR), and stellar 
mass. We find that there is a smooth transition in the clustering of galaxies as a function of this color 
from weak clustering among blue galaxies to stronger clustering for red galaxies. The clustering of 
galaxies within the "green valley" has an intermediate strength, and is consistent with that expected 
from galaxy groups. The results are robust to the correction for dust extinction. The comparison 
with simple analytical modeling suggests that the halo occupation number increases with older star 
formation epochs. When splitting according to SSFR, we find that the SSFR is a more sensitive tracer 
of environment than stellar mass. 
Subject headings: ultraviolet: galaxies 



1. INTRODUCTION 

The nature of star forming galaxies has deeply changed 
with the evolution of the Universe. Besides the strong de- 
crease in the cosmic star formation rate since z ~ 1 (e.g., 
iHopkins fc Beacoml[200llSchiminovich et al.ll2005h . the 
bulk of star formation has shifted from high stel- 
lar mass to low stellar mass systems, ( "downsizing", 
see jBundv et all 120061 : ICowie et alJl!996t Uuneau et all 
l2005h , a phenomenon w hich exhibits a si milar trend in 
dark matter halo mass (|Heinis et al.l l2007h . In the lo- 
cal Universe results from optical surveys show that the 
galaxy population is bimodal, divided at a stellar mass 
of ~ 3 x 1O 1O M0 between actively star-forming blue 
sequence galaxies and passively evolving red sequence 
galaxies (|Kauffmann et aLll2003al ). This bimodality in 
galaxy propertie s has been obse r ved from low redshift 
up to z ~ 15 dBell et al.l 12004 ICirasuolo et all 120071: 
ICassata et all 120081; iFranzetti et all 120071 ) although the 
proportion between the red and blue sequence has been 
changing. Several studies indeed show that the den- 
sity of red galaxies has been increasing since z = 2 
while the density of galaxies on the blue sequences re- 
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main ed roughly constant (|Bell et alj 1200 4: Fabe r et al.l 
2007). These are further pieces of evidence that some 
previously star-forming galaxies have quenched their star 
formation and moved from the blue to the red sequence 
at a rate balanced by ongoing star formation among the 
remaining blue se quence galaxies (jArnouts et all [20071 : 
iMartin et al. 2007). Various mechanisms can be involved 
in the quenching of star for mation. Galaxy me rgers 
can trigger AGN activity (e.g.. lHopkins et al]|2007f ). and 
hence AGN feedbac k, which is efficien t in further pre- 
venting gas cooling (jCroton et al.ll2006l ). Supernovae ex- 
plosions inject energy into the interstellar medium and 
may expel gas fro m galaxies in low- mass halos (e.g., 
iDekel &: Silld fT9 86L Galaxy environment is also an im- 
portant factor in this evolution: in clusters for instance , 
several proce sses such as starvat ion (jLarson et alj [l980). 
haras sment (|Moore etaLl Hflflflj) or ram pressure strip- 
ping (|Gunn fc GottJll972l ). which can be effective a t large 
distances from the cluster center ( Ch ung et al.l 120071 : 
iTonnesen et al.ll2007h . may combine to lead to morpho- 
logical transformations and removal of the gas reservoir. 
Numerical studies also suggest that accounting for grav- 
itational heating through satellite accretion can repro- 
duce the ob served trends between star formation and en- 
vironment (|Khochfar fc Ostriker|[2008l ). 

A common tool used to place constraints on the phys- 
ical mechanisms involved in this evolution has been the 
study of galaxies' properties and their link with environ- 
ment within the color-magnitude diagram. Recent stud- 
ies have showed that at fixed luminosity the fraction of 
red galaxies increases in denser environments, while the 
mean color of the red a nd blue sequences ar e fairly insen- 
sitive to local density ( jBalogh et al.ll2004T ). iHogg et ail 
(2003) found that for blue galaxies, at fixed star forma- 
tion history, stellar mass is weakly correlated with envi- 
ronment, and that lumino us and faint red galax ies reside 
in denser environments. ICucciati et al.l (|2006l ) showed 
that the relation between color and local density evolved 
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between z = 1.5 and z = 0, which suggests that the links 
between galaxy properties and environment observed at 
the present epoch are not only determined by galaxy for- 
mation conditions ("nature") but also by contributions 
from galaxy evolution processes ( "nurture" \ 

Ex tendin g these studie s usin g GALEX ([Martin et al.l 
12005ft data. IWvder et al.l (|2007l ) showed that the galaxy 
distribution is less bimodal in ultraviolet (UV) optical 
color-magnitude diagrams, owing to the greater sensi- 
tivity of the UV to star formation and dust. Previous 
work also showed that there is a very tight correlation 
between the (NUV — r) a b s 10 color and star formation 
history (jSalim et al.ll2005h . UV-optical colors further- 
more highlight a galaxy population within the interme- 
diate region between the red and the blue se quences, the 
s o-called "green val l ey" (TMartin et all 12007ft . 

iSchawinski et al.1 (|2007ft studied early-type galaxies 
within the UV-optical color magnitude diagram; they 
showed that this population exhibits residual star for- 
mation, and that the fraction of UV-bright early-type 
galaxies is higher in lower density environments. In 
this paper, we extend this study by using the combina- 
tion of GALEX imaging with SDSS spectroscopic data 
to measure the dependence of spatial clustering on the 
(NUV — r) a bs color, in order to link star formation his- 
tory with environment, the latter being traced by the 
large-scale clustering. We use here a simple power-law 
approximation to describe the clustering of the galaxies, 
in order to estimate the mass of the host dark matter 
halos. 

The outline of this paper is as follows: in Section [2] we 
present the samples from the GALEX surveys, as well as 
the internal dust correction procedure. Section [3] deals 
with the clustering analysis tools, and contains a discus- 
sion of the systematic errors previously n oticed for clus- 
terin g measurements from GALEX data (Milliard et al. 
12007t h We present the results in Section [4] discuss them 
in Section and conclude in Section O Throughout this 
paper we use the concordance ACDM cosmology with 
f2 m = 0.3 and J1a = 0.7; we note Hq — 100/ikm s _1 
Mpc -1 . We use h = 1 for clustering computations to fa- 
cilitate comparisons with other studies, and h = 0.7 for 
all absolute magnitudes derivations. 

2. DATA SAMPLES 

We consider observations from two different GALEX 
surveys: the Medium Imaging Survey (MIS) to use the 
UV photometry for the overall galaxy population de- 
tected by the SDSS spectroscopic survey, and the All sky 
Imaging Survey (AIS ) , to build a UV selected sample. 
IBudavari et al.l (|2009f ) described in detail how we built 
the catalogs and footprints; we only summarize here the 
main characteristics. 

We computed the unique area observed in the GALEX 
MIS and AIS. The fields target pre-defined positions on 
the sky. These 47,162 points define disjoint cells by in- 
cluding in each cell those regions of sky closer to the cen- 
ter of that cell than to the center of any other cell, similar 
to the way in which a Voronoi tessellation is defined. The 
intersection of the circular field-of-view of a given field 
and the targeted cell is the primary region for this field. 
We use here only objects within the primary resolution 

10 We denote absolute magnitudes by the 'abs' subscript. 



of AIS or MIS fields. The union of the primary region 
from all fields defines the pri mary footprint of a surv ey. 
We then use the formalism of IBudavari et al.1 (|2008ft to 
compute the footprint of the GALEX-S DSS overlap. 

Th e GALEX 3 rd Release (GR3; iMorrissev et~all 
12007ft dataset was cross-matched to the SDSS DR6 
(|Adel man-M cCart hv~et al.l [2008) using a matching ra- 
dius of 4". We used the GALEX sources with one 
and only one SDSS counterpart within the search ra- 
dius. In the following we only consider SDSS primary 
sources. We further restricted the SDSS sources to 
spectroscopically confirmed galaxies with z < 0.3 and 



14.5 < 



Tpetro 



< 17.6. The bright cut avoids the shred- 



ded galaxies, while the faint selection ensures a homo- 
geneous spectroscopic completeness over the sky. For 
GALEX sources, we retained the sources within 0.5 de- 
grees of the field centers, as the photometry accuracy 
decreases beyond this limit, and artifact s concentrate at 
the e dges of the GALEX field of view (|Morrissev et al.1 
[2005ft . We censor sources using GALEX (flagmap values 
2 and 4) and SDSS (BLEEDING, BRIGHTJ3TAR, TRAIL and 
HOLE) masks. We use the intersec tion of the GALEX an d 
SDSS footprints as described in IBudavari et alJ |2008). 
We corrected GALEX apparent magnitudes for Galac- 
tic extinction using the lSchlegel et ail (1 19981 ) dust maps 
and the formulas described in IWvder et alJ (|2007t ) . We 
estimated fc— corrections 11 at z — after fitting the 
7- band photometry (FUV, -NUV, u,g,r 7 i,z) to a set 
of iBruzual fc Chariot) ([2003ft templates 12 ; the absolute 
magnitudes are fc— corrected to z — unless otherwise 
stated. We also derived fc— corrections at z = 0.1 by 
shifting the filters blueward in wavelength by a factor 
1/(1 + 0) = 1/1.1, in order to compute dust extinction 
correction (see 



2.1. MIS/SDSS selection 

The MIS observations are deep enough to probe the 
overall popula tion of the SPSS M ain Galaxy spectro- 
scopic sample ([Strauss et al.ll200"2ft . We used MIS fields 
with exposure times greater than 1000s for the NUV 
band, and objects with NUV <23. Fig.[T]shows the color- 
magnitude diagram in the (NUV, NUV-r) plane for the 
MIS sources. The dashed line corresponds to the faint 
cut in r, and the dotted line to the bright cut we adopt. 
This sample is obviously not UV-selected, as the density 
decreases for objects fainter than r = 17.75, the SDSS 
target magnitude limit. There is however a clear hint of 
a bimodality in color, with a separation around NUV-r 
~ 4.75 between the red and the blue sequence. The final 
sample contains 21,895 galaxies over 598.1 sqdeg. 

2.2. AIS: UV selection 

We used AIS observations to build UV-selected sam- 
ples from the GALEX-SDSS spectroscopic samples. Fig- 
ure [2] shows as before the color-magnitude diagram for 
AIS objects, with an additional dotted line representing 
the apparent magnitude cut at NUV=19 we use to ensure 
a UV-selected sample. The fraction of objects brighter 
than NUV — 19 and fainter than r = 17.6 that are 
lost is 3%. We note that the NUV cut removes most of 
the galaxies on the red sequence; hence our UV-selected 
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Fig. 1. — Color magnitude diagram for MIS selection. The 
dotted and dashed lines show the cuts in r adopted for SDSS 
spectroscopic objects. 

galaxy sample mainly contains blue galaxies. The size 
of the sample is 14,351 galaxies, over an area of 3645.8 
sqclcg. 
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Fig. 2. — Color magnitude diagram for AIS selection. The 
dotted and dashed lines show the cuts in r adopted for SDSS 
spectroscopic objects, and the solid line the NUV cut chosen 
to ensure an UV-selected sample. 



2.3. Dust extinction correction 

In the following we study the dependence of clustering 
as a function of the {NUV — r) a b s color, which correlates 
tightly with star formation history and provides a very 
good separation between galaxy populations. However 
interstellar dust absorbs and scatters stellar light, red- 
dening galaxy colors and interfering with the interpreta- 
tion of the colors in terms of the star formation history. 
Several recipes are available to correct for dust redden- 
ing, but one often has to rely on various assumptions 
about the dust content of galaxies. For instance, galax- 



ies with intermediate colors may have redder colors than 
blue galaxies because of dust, or because of older stellar 
populations. Hence applying a given dust extinction cor- 
rection may imply potential errors for transition or red 
galaxies. We will therefore present clustering measure- 
ments as a function of both non-correcte d and corrected 
colors . We use the relation presented by iJohnson et al.l 
(200<|: 



A FUV = 1.27 - 1.56x + 1.35?/ - 1.2Axy 
ir = D n (4000) - 1.25 
y={NUV-r) absfi . 1 -2 



(1) 



to esti mate the extinc t ion c orrection in the FUV band, 
Afuv- IJohnson et alJ (|2006h used t he nar rowband color 
£> n (4000) from iKauffmann et all (|2003af ) to calibrate 
their relation, while we consider here the Z?(4000) index 
directly available from the SDSS DR6 pipeline, which 
involves the ratio of the flux over a longer wavelength 
range. We checked that for objects in common between 
the SD SS DR6 and the sample used bv[Kauffmann "et al.l 
(|2003al ) that the difference between their estimate of 
Z?„(4000) and our estimate of Z?(4000) is reasonably 
small with a dispersion of ~ 0.06. Finally, we did not 
apply any dust correction to objects without star forma- 
tion, defined here as those galaxies with EW(Ha) < 0. 
These objects represent 7% of the MIS/SDSS sample and 
less than 1% of the AIS UV-selected one. Negative values 
of Afuv accor ding to eq. Q] ar e set t o 0. 

We used the ICalzetti et alJ (|2000| ) extinction law to 
convert the extinction in the FUV band into extinction 
in the NUV, g and r bands with: Anjjv = OSlApuv, 
fuv an( l A r = 0.35 Afuv- 



A g = 0.46A 



3. CLUSTERING ANALYSIS 



We chose to quantify galaxy clustering in our samples 
using the real-space correlation function. As redshift sur- 
veys suffer from redshift distortions due to peculiar ve- 
lociti es, we followed the common practice ([Fisher et al.1 
fl99l of measuring first the redshift-space correlation 
function £(rp,7r), which is useful for distinguishing be- 
tween the contribution of galaxy pairs along the line of 
sight 7r and the dire ction perpendicu l ar to the line of 
sight r p . We used the lLandv fc S zalav (1993) estimator: 



£(?>,7r) 



DD - 2DR + RR 
RR 



(2) 



where DD, DR and RR are the data-data, data-random 
and random-random pair numbers, respectively. To gen- 
erate random catalogs, we distributed objects randomly 
on the sky within the GALEX-SDSS footprint in non- 
masked areas, which takes care of the angular selection 
function. We then randomly picked galaxy properties 
(such as redshift and magnitude) that we assigned to the 
random objects. This procedure has been shown to pro- 
vide accurate results, and enabled us to study the depen- 
dence of clustering on galax y properties by fully taking 
into account selection effects ()Li et al.ll2006al ). While the 
GALEX-SDSS overlap is fairly patchy, we checked that 
the integral constraint is negligible for this study. 

We used 50 times more random objects than the num- 
ber of galaxies in the sample under consideration to com- 
pute the correlation function. We found out that this is 
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the minimum factor to estimate accurately the correla- 
tion function on small scales, given the fairly low statis- 
tics of the samples we use here. 

As redshift distortions occur only on the line-of-sight, 
one can get rid of their effects by computing the projected 
correlation function 



w p (r p ) = 2 



(3) 



After checking for convergence of the integration step, 
we used here ir max = 25ft _1 Mpc. The projected cor- 
relation function is finally related to t he real-space cor- 
relati on function £(r) by the relation (|Davis fc Peebles] 
1983) 



w p {r p ) = 2 £ (Jrl + y 2 ) dy 



(4) 



where y is the real-space distance along the line of sight. 

Under the assumption of a power-law £(r) = (V/ro) 7 
this relation yields 



w p {r p ) 



r(*)r(V) 
r(3) 



(5) 



We will use the parameters r*o, the correlation length, 
and 7, the slope of the real-space correlation function, 
to describe the clustering of the galaxies. We use jack- 
knife sampling to derive the error bars on w p (r p ); this 
procedure has been s hown to provide acc urate error mea- 
surements on w p (r p ) ([Zehavi et alj|2005l) . Specifically we 
divided our samples into N separate regions (32 for AIS 
sample, 20 for MIS), large enough so that their volume is 
significant. We then computed w p (r p ) for the full sample 
but removing one jackknife region in turn. For a given 
bin r p -, the errors from these N samples are then given 
by: 



N-l N 



3=1 



(6) 



3.1. Fiber collision correction 



Due to the finite size of the fibers used to take the SDSS 
spectra, fibers can not be set closer than 55", resulting 
in a lack of close pairs at comoving scales smaller than 
~ 80/1" 1 kpc (at z = 0.1). We use a statistica l appr oach 
to correct from this, following lHawkins et alj (|2003f) and 
ILi et all (l2006ah . 

We measured the angular correlation function w t (0) for 
the objects classified as galaxies and selected for spectro- 
scopic measurement by the SDSS (i.e. target galaxies, 
open circles on fig. [3£i) and for SDSS objects with as- 
signed galaxy spectra w z (9) (open squares on fig. [3^,) 
within the AIS footprint (large enough to provide an ac- 
curate measurement). We then used the ratio 



F{6) 



l + w t (0) 
l + w z {6) 



(7) 



as an estimate for the correction factor needed to correct 
for fiber collision. The measured values of this ratio are 
shown in fig. [3Jd); we used a functional form for this cor- 
rection factor represented by the solid lines on fig. [HJd. 
F(9) can be then used to weight the galaxy pairs accord- 
ing to their angular distance. As a test we recomputed 



the angular correlation function of the SDSS objects with 
assigned galaxy spectra using this weight (filled triangles 
on fig. [3^,); the agreement between this result and the 
angular correlation function of the target galaxies is ex- 
cellent. 

3.2. Systematic errors 

In iMilliard eTall (|200l we noticed that the angular 
correlation function from GALEX data shows signatures 
similar to those expected in the case of low-level field-to- 
field variations. This prevented us from estimating the 
angular correlation function from the whole dataset at 
once. The approach we used was to estimate the cluster- 
ing from the angular correlation function by measuring 
w(8) independently for each GALEX field, and then av- 
eraging the results. This method cannot be used in the 
specific case of the GALEX-SDSS spectroscopic samples, 
given the low statistics. Hence we performed a series of 
tests to see whether the order of magnitude of the effects 
on the spatial correlation function are the same as those 
on the angular correlation function. As in Milliard et al. 
(|2007f ). we use the Mann- Whitney statistic to determine 
which GALEX fields should be considered as outliers in 
terms of UV photometry. We consider potential out- 
liers fields those with a p- value from the Mann- Whitney 
test lower than 0.05 13 . We show fig. [4] (left) the pro- 
jected correlation functions w p {r p ) (computed using all 
data at once) of the full AIS sample 14 (filled circles), of 
the fields expected to be outliers (p-value < 0.05, filled 
squares, 16% of the full sample) and the fields that are 
expected to have an homogeneous UV photometry (p- 
valuc > 0.05, filled triangles). The correlation function 
of the outlier fields is very different from the fields with a 
p-value > 0.05, showing a higher amplitude, and a shal- 
lower slope. Fitting these results by power laws (see % 2 
contours on fig. |4j right) shows that the outlier fields 
have a correlation length roughly 1.5 times larger than 
the fields with high p-values. These outlier fields have 
a fairly small impact on the clustering parameters de- 
rived for the full sample; excluding or including them 
from the analysis does not change the conclusions of this 
study. Nevertheless, these fields indeed show a very dif- 
ferent clustering. This stronger clustering may be spuri- 
ous, or may be real. In the latter case, it might be due 
to large-scale structure. To further test this hypothesis, 
we computed w p (r p ) from the entire SDSS data within 
the AIS footprint, splitting the fields as before, based on 
the Mann- Whitney results obtained from UV photome- 
try. The trends are similar, with the SDSS data within 
the footprint of the AIS outlier fields showing a signifi- 
cant stronger clustering 15 . We also performed the Mann- 
Whitney test on SDSS data, using the u band; 60% of 
the outlier fields are the same as the ones according to 
NUV photometry. The clustering derived from the fields 
classified as outliers with both u and NUV photometry 
actually dominate the clustering signal observed for the 
outliers on fig. |H 

13 This means that the probability that the photometry of these 
fields comes from the same parent distribution as the other fields 
is lower than 0.05. 

14 The trends are similar for other apparent magnitude cuts and 
the MIS sample. 

15 Trends are also similar for SDSS data within MIS footprint, 
when splitting fields according to MIS photometry. 
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Fig. 3. — Fiber collision correction. Left: angular correlation function of SDSS target galaxies (wt(0), open circles), SDSS 
sources with assigned galaxy spectra (w z (9), open squares). The filled triangles show the angular correlation function of the 
SDSS sources with assigned galaxy spectra computed using the correction from fiber collision. Right: the symbols show the ratio 
F(9) = [f + w t (6)]/[l + w z (9)] derived from the measurements shown on the left plot. The solid line shows the functional form 
we used for this ratio in order to correct for fiber collisions. The dotted line shows the 55"limit. 

As a final test, a friend-of-friend structure finding al- 
gorithm applied on SDSS spectroscopic data shows that 
there is a slight trend for fields with lower p-values de- 
rived from UV photometry to have larger galaxy groups, 
as recovered by the friend-of-friend algorithm. These 
tests indicate that the stronger clustering, as measured 
by the combination of GALEX and SDSS spectroscopic 
data, observed for the fields considered as outliers based 
on the Mann- Whitney test applied on UV photometry 
is genuine. We checked for instance that the large-scale 
struc ture known as the "Sloan Great Wall" (Gott ct al. 
I2005D is not within our sample; our tests rather show that 
the fields we classified as outliers contain galaxy groups 
with larger numbers of members. In consequence, these 
fields need to be considered in the global analysis. We 
hereafter use the full AIS and MIS samples, and com- 
pute the projected correlation function from the whole 
dataset. 



4. RESULTS 

4.1. Redshift space clustering 

Figure [5] shows the redshift space correlation function 
£(r p ,7r) computed from the MIS sample (left) and the 
AIS sample (right). The contour s show the w ell-known 
effects due to peculiar velocities (|Kaiserl fl987): at small 
scales, galaxies in clusters lead to distortions along the 
line of sight (fingers of God); this effect is clearly no- 
ticed in our MIS/SDSS sample. As expected, this feature 
is much less prominent in the UV-selected AIS sample, 
which contains very few red galaxies. At large scales, 
coherent infall on large structures causes a squashing of 
the contours, clearly seen in both samples. 

4.2. Clustering as a function of (NUV — r) a b s color 

We first consider the MIS sample in order to use the 
very good separation between red and blue sequence 
galaxies in (NUV — r) a b s color. The distribution of the 
galaxies in the UV-optical color diagram is shown fig. 



6l As noticed by previous studies (|Martin et al.1 120071 : 
Wvder et al.ll2007l ). the color distribution shows a clear 
bimodality between the blue and red clouds and a tran- 
sition population located in the so-called green valley. 
The distribution of the galaxies in the UV-optical color 
diagram after applying the dust extinction correction is 
shown in fig. The shape of the distribution is quite 
different once this correction is applied, with a less promi- 
nent red sequence. 

We divided our sample into four bins of (NUV — r) a b s 
color in the following way. We set up first the green valley 
for both versions of the color magnitude diagram (3.75 < 
(NUV— r) a bs < 4.5 without extinction correction; 2.75 < 
(NUV— r) a bs < 3.5 with extinction c orrection), followin g 
previous studies (see fig. 9 and 25 of I Wvder et ai1l2007f ). 
This also defines a natural cut for the red and the blue 
populations. We further divided the blue cloud by the 
median color of this population. The cuts we adopted 
are shown as dotted lines on figures [5] and \7\ 

The correlation functions are shown fig. [8] (for color 
cuts without dust extinction correction) and fig. [9] (for 
color cuts with dust extinction correction) , as well as the 
power-law fits. In both cases, there is a monotonic in- 
crease of the clustering amplitude for redder galaxies, as 
well as a steepening of the correlation function. The re- 
sults obtained with (see clustering parameters table \2§ 
or without the dust extinction (table QJ correction are 
very similar, while the only criterion to define the color 
cuts has been to isolate the green valley. This is ex- 
pected to some extent as the extinction correction we ap- 
plied is a function of the uncorrected (NUV — r) a b s color 
and D(4000), which are themselves fairly well correlated 
with one another, (e.g. iMartin et al.ll2007t IWvder et al.1 
I2007D . Note also that the dust corrected green valley 
overlaps with the non dust corrected one at the 50% 
level. One could argue that the intermediate clustering 
of the green valley galaxies we select with non-dust cor- 
rected colors is due to a mix of dusty star forming galax- 
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r p [h 1 Mpc] r [h 1 Mpc] 

Fig. 4. — Left: projected correlation function w p (r p ) from the AIS sample with NUV< 19. Circles show w p (r p ) for the 
full sample (All), squares w p (r p ) for the fields flagged as outliers in terms of UV photometry by the Mann- Whitney statistics 
(p < 0.05), and triangles w p (r p ) for the remaining fields (p > 0.05). Solid lines represent the best power law fits. Right: \ 2 
contours of the power law fits in the (ro, 7) plane; the contours labels include the legend for each correlation function of the left 
panel plot. Dashed lines show the contour at xLn + 1; while the thin and thick solid lines show the contours at 68% and 90% 
respectively. 




Fig. 5. — Redshift space correlation functions. Left: Contours of £(r p ,7r) for the MIS sample. Right: Same as left panel but for 
the AIS sample. The solid black line shows the contour at £(r p ,7r) = 1. The contour levels are logarithmically spaced. 



ies and galaxies with older stellar populations. However, 
the fact that this result holds when using dust corrected 
colors suggests that dust does not induce such a large 
population mix in this region of the color-magnitude di- 
agram. Another possibility is that the gr een valley is 
actual ly an intrinsic composite population; ISalim et al.l 
(2007) showed that it is built from normal star forming 
galaxies, star forming/AGN composite, and pure AGNs. 
However, tests show that within the green valley the 
large-scale clustering of AGNs and non AGNs is very 
similar (see sec. 15. 1|) . 

We also used the AIS sample to investigate the clus- 
tering as a function of the (NUV — r) a b s color within an 



UV-selected sample. The color distribution for this sam- 
ple is shown on the right panels of figures [5] and [7] with 
dashed lines. We used the same color cuts as for the 
MIS/SDSS sample and also define three cuts in order to 
get similar subsample sizes; we use here colors corrected 
from dust extinction. The correlation functions are pre- 
sented in fig. [TU] and the clustering parameters in table 
[3J As observed in our MIS/SDSS sample, the correlation 
length increases for redder galaxies, although our UV- 
selected sample is built mainly from galaxies on the blue 
sequence (see fc|2.2|) . This increase is not observed if we 
use the same cuts as for the MIS/SDSS sample. The AIS 
correlation function for the 1.85 < (NUV — r) a b s < 2.75 
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TABLE 1 

MIS/SDSS SAMPLES CLUSTERING RESULTS* 



Sample N gal (z) r [h 1 Mpc] -y b 8 

-1.0 < (NUV - r) abs < 2.8 6293 0.08 2.93 ± 0.09 1.58 ± 0.05 0.66 ± 0.02 

2.8 < (NUV - r) abs < 3.75 6291 0.10 4.0 ±0.12 1.67 ± 0.03 0.86 ± 0.02 

3.75 < (NUV - r) abs < 4.5 2587 0.10 5.48 ± 0.21 1.77 ± 0.04 1.13 ± 0.04 

4.5 < (NUV - r) abs < 8.0 6702 0.10 6.93 ± 0.32 1.84 ± 0.03 1.43 ± 0.06 

Without dust extinction correction 



TABLE 2 

MIS/SDSS SAMPLES CLUSTERING RESULTS* 



Sample 






<*} 


r [h' 1 Mpc] 


7 


bs 


0. < (NUV - r) abs , corr 


< 1.85 


6635 


0.08 


3.0 ± 0.1 


1.60 ± 0.05 


0.67 ± 0.02 


1.85 < (NUV - r) abs , corr 


< 2.75 


6534 


0.10 


4.25 ± 0.12 


1.67 ± 0.03 


0.89 ± 0.02 


2.75 < (NUV - r) abs , corr 


< 3.5 


2769 


0.10 


5.25 ± 0.24 


1.80 ± 0.05 


1.09 ± 0.05 


3.5 < (NUV - r) abs , corr 


< 7.0 


5816 


0.10 


7.07 ± 0.32 


1.87 ± 0.03 


1.46 ± 0.06 



With dust extinction correction 
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Fig. 6. — Left: color-magnitude diagram from the 
MIS/SDSS sample. The dotted lines show the cuts used 
to measure the projected correlation function. Right: his- 
tograms of the (NUV — r) abs color, with the same cuts shown 
on the left panel; solid line shows MIS/SDSS sample, and 
dashed line AIS UV-selected sample. 

bin is noisy because of small statistics. 

Interestingly, our results are poorly fitted by a simple 
power-law. In particular at scales r p < O.l/i -1 Mpc, the 
correlation functions for the blue galaxies in both the 
MIS and AIS samples show an excess of pairs, while on 
the contrary the correlation function of our reddest bin 
shows a lack of satellites (figs. [8] and [9]). We note that 
for the red galaxies, we are biased against low masses 
and low specific star formation rates objects (see fig-fTS]). 
These faint red galaxies a re known to exhibi t a strong 
clustering at small scales ([Zehavi et al.l 12005). For the 
blue galaxies, our results suggest a higher satellite frac- 
tion for bluer UV galaxies, which is similar to what has 




Fig. 7. — Same as fig. [6] but using dust extinction corrected 
absolute magnitudes. 

been observed a t high redshift from UV-selected samples 
(|Lee et al.ll2006l: lOuchi et al.ll2005h . 

We show in fig. [TT] the combination of our results for 
both the AIS (filled and open squares) and MIS samples 
(filled circles). Thanks to the cut in (NUV — r) a b s color, 
we are able to probe a large dynamic range in terms of 
clustering strengths, with values of the correlation length 
ro increasing from ~ 1.8/i~ 1 Mpc for the bluest galaxies, 
to ~ 7.1ft. -1 Mpc for the reddest. The UV-selected and 
MIS/SDSS samples follow the same trends. However the 
purely UV-selected galaxies from the AIS survey show a 
weaker clustering strength at a given color. The corre- 
lation length from the AIS sample is lower than the one 
obtained from the MIS/SDSS sample if the same cuts 
are applied. The mean color for the AIS subsample with 
< (NUV — r) a bs < 1-85 is lower than the correspond- 
ing one for the MIS, but the conclusion holds if different 
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TABLE 3 

UV-SELECTED SAMPLES CLUSTERING RESULTS 



Sample N gal (z) r [h~ x Mpc] 7 6 8 

0. < {NUV - r) abSjCOII < 1-85 10861 0.05 2.59±0.11 1.69 ± 0.05 0.57 ± 0.02 

1.85 < {NUV - r) abs ,corr < 2.75 2596 0.06 2.U+g_% 1.51 ± 0.12 0.56±q 07 

0. < (NUV - r) abS:Corr < 1.35 4507 0.05 1.79 ± 0.17 1.55 ± 0.10 0.45 ± 0.04 

1.35 < (NUV - r) abs corr < 1.64 4503 0.05 2.07 ± 0.22 1.51 ± 0.10 0.51 ± 0.04 

1.64 < {NUV - r) abs COII < 4. 4592 0.05 2.75 ± 0.16 1.54 ± 0.06 0.62 ± 0.03 



1000 




1 I 

0.01 



4.5 <(NUV-r) abs <8. * 
3.75<(NUV-r) abs <4.5 * 

2.8 <(NUV-r) abs <3.75 ■ 
-1. <(NUV-r) abs <2.8 • 
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Fig. 8. — Projected correlation functions of MIS/SDSS sam- 
ples cut by {NUV — r) abs color (without dust extinction cor- 
rection). 



1000 : 



100 : 



10 r 



1 L_ 
0.01 




3.5 <(NUV-r) abs , 
2.75<(NUV-r) abs-corr <3.5 * 
1.85<(NUV-r) abscorr <2.75 ■ 
0. <(NUV-r) abs , OT ;<1.85. 
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Fig. 9. — Projected correlation functions of MIS/SDSS sam- 
ple cut by {NUV — r) abs color (with dust extinction correc- 
tion). 

cuts are applied and the mean color is similar. 

4.3. Comparison with previous studies 

To assess the impact of the NUV selection on our 
MIS/SDSS sample, we can directly compare our re- 
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Fig. 10. — Projected correlation functions of AIS UV- 
selected samples in bins of {NUV — r) abs color. 



suits with those of IZehavi et al.l (|2002), who divided 
their sample between blue and red galaxies using the 
limit {u — r) a bs — ^-8- They found for red galaxies ro = 
6.78±0.23h~ 1 Mpc and 7 = 1.86±0.03; for blue galaxies 
r = 4.02 ± O^/i^Mpc, 7 = 1.41 ± 0.04. Applying the 
same cuts to our MIS/SDSS sample, we obtain for red 
galaxies r Q = 5.20 ± 0.15/1" 1 Mpc, 7 = 1.77 ± 0.02; for 
blue galaxies r = 3.24 ± 0.17 /t" 1 Mpc, 7 = 1.66 ± 0.07. 
We find the same trends as IZehavi et al.l ((2002) with 
lower correlation lengths. 

We observe a steepening of the slope from 7 = 1.58 ± 
0.05 for the bluest galaxies to 7 = 1.84 ± 0.03 for the 
reddest for the results from the MIS/SDSS sample with- 
out the dust extinction correction. This increase of the 
slope with color is similar to what has been already ob- 
served by previous studies from optical samples at low 
rcdshift. Selecting active and p assive galaxies from their 
spectra, iMadgwick et alj (|2003f ) found that they exhibit 
power laws with respective slopes 7 = 1.60 ± 0.04 and 
7 = 1.95 ± 0.03. iBudavari etai] (|2003h observed simi- 
lar trends in their angular correlation function measure- 
ments as a function of spectral type. Their blue samples 
have a mean slope 7 = 1.67 ±0.11 and their red samples 
have 7 = 1.99 ± 0.06. 

The correlation length we obtained from the 
MIS/SDSS sample for the bluest galaxies in {NUV - 
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Fig. 11. — Dependence of the correlation length ro on the 
(NUV — r) a b s ,corr color. The error bars on color are the stan- 
dard deviation of the colors in each bin. Filled circles show 
the MIS/SDSS sample results. Squares show the AIS (UV- 
selected sample) results: open squares using the same cuts as 
for the MIS/SDSS sample, and filled squares cuts defined in 
order to have similar number of galaxies per subsample. 

Oabs °f r o ~ 3.0 ± 0.1/i _1 Mpc is lower than the values of 
3.6 — 4.5/i _1 Mpc foun d for blue galaxies in optical sam- 
ples s uch as the SDSS (Budavari ct al. 2003; Zeha vi et 
120051) or 2dF (|Madgwick et a l. 2003; Norberg et aT1l2002ft 
at similar redshifts. This lower value is due to the com- 
bination of our selection which biases against low stellar 
mass and low specific star formation rate objects, and 
the larger star formation history dynamic range probed 
by the (NUV — r) a b s color as well. 

The ro value we measured for the bluest bin of our 
UV-selected sample, r = 1.79 ± O.m^Mpc, is among 
the lowest observed in the local Universe. It is slightly 
lower than the value r = 2.31± l ;!/i" 1 Mpc we obtained 
from a flux-limited (NUV<21.) UV-selected sample 
(|Milliard et al.ll200l . 

These estimates are closer to the value s derived for 
gas-ric h, H-I selected galaxies obtained by iMever et al.l 
pOOl of r = 2.70 ± O^lfr^Mpc and 7 = 1.56 ± 0.1. 
This correspondence might indicate that we can asso- 
ciate the bluest galaxies with gas-rich galaxies, which is 
in agreement with the fact that the gas fraction o f blue 
galaxies is larger (see, e.g.. ISchiminovich et al.ll2007l ). 

Transition galaxies in the green valley indeed exhibit 
an intermediate clustering pattern, although comparison 
with other studies might lead to a different conclusion 
depending on the selections used. For instance the cor- 
relation leng t h we measure is close to that derived by 
iZehavi et al.l (|2005t ) for red galaxies selected by g — r 
color and -20 < r abs < -19 (r = 5.7/i _1 Mpc). On the 
other hand, our ro value for the green valley is also sim- 
ila r to one of the f aint s amples of blue galaxies studied 
bv lBudavari et alJ (|2003l ) (T3 sample, r = 5.9±0.51h- 1 
Mpc), who selected galaxies acc ording to their sp ectral 
energy distribution. At z ~ 0.85 iCoil et al.l (|2008l) mea- 
sured the projected correlation function of green galaxies, 
selected using the (U — B) restframe color. The correla- 
tion length of ro = 5.17 ± 0.42/i~ 1 Mpc that they derived 




D(4000) 

Fig. 12. — Extinction correction in the FUV band, Afuv, 
as a function of the 4000 A break, D(4000) for the MIS/SDSS 
sample. The extinction correction is obtained from eq.[T] The 
three color shades code objects on the blue sequence, green 
valley and red sequence, from left to right. 

is similar to the one we obtained. 

5. DISCUSSION 

5.1. Intermediate clustering strength of green valley 
galaxies 

The results presented in sect. |4] show that green valley 
galaxies have an intermediate clustering strength com- 
pared to bluer or redder galaxies. This holds whether 
dust extinction is applied or not. As the green valley 
is a mix of intermediate star formation history galax- 
ies and dusty star forming ones, interpretation of this 
result relies on the efficiency of the dust correction pro- 
cedure to distinguish between these populations. Fig- 
ure [T2] shows the extinction correction in the FUV band, 
Afuv as a 

function of the 4000 A break, D(4000) for 
the MIS/SDSS sample. Galaxies within the green valley 
(3.75 < (NUV - r) abs < 4.5) lie in the middle (green) 
region of this plot. This figure shows that the dust cor- 
rection prescription we used takes into account the evo- 
lutionary stage of the galaxy population. For instance, 
more evolved galaxies with high D(4000) or low specific 
star formation rates have a lower extinction correction 
than younger (low D4000) or more active galaxies (high 
specific star f ormation rates ) , in a greement with recent 
analysis from ICortese et all (|2008f ) . While still subject 
to some uncertainties, this procedure should clean up 
the green valley population from the highly dusty popu- 
l ation. 

iMartin et alJ (|2007l ) noted that the fraction of galaxies 
with AGNs (in their case type 2, or narrow-line AGNs) 
peaks in the green valley. There is no study devoted to 
the clustering of type 2 AGNs to whic h we can directly 
compare our results. iLi et alJ (|2006bl ) studied however 
the cross-correlation of narrow-line AGNs with galaxies 
within SDSS. In particular, they showed that their clus- 
tering is in agreement with the assumption that a high 
fraction of them are within central galaxies in dark mat- 
ter halos, and that they reside mainly within halos with 
masses 10 12 -10 13 M Q , values consistent with our results 
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Fig. 13. — Projected correlation functions for MIS green val- 
ley galaxies, assuming dust extinction correction. The sample 
is here restricted to the overlap with SDSS DR4, and galax- 
ies divided between galaxies with (filled circles) or without 
AGNs (filled squares). 



for the green valley (see ^5 . 2[) . Note also that several 
studies pointed out the "intermediate" nature of galaxies 
hostin g AGNs (|Heckman et al.1 l2004t iKauffmann et all 
2003b). For instance, the average spectra of AGNs galax- 
ies host is remarkabl y similar to those of disk galaxies 
with Sb Hubble type (|Yip et al.ll2004al|b1 ). 

The intermediate clustering we measured for transi- 
tion galaxies may be then mainly driven by galaxies 
hosting AGNs. To test this further, we restricted our 
MIS/SDSS sample to the overlap with SDSS DR4, in or- 
der to discriminate, within the green valley, AGNs and 
non AGNs, based on the catalogs from IKauffmann et all 
(j2003cj ). IKauffmann et all (|2003ch selected type 2 AGNs 
from star form i ng gal axies using the diagram proposed by 
iBaldwin et al.l (I1981D. We note that th e objects classified 
as 'AGNs' bv lKauffmann et al.l (|2003d ) include star form- 
ing/ AGNs composites; the low numbers of this sample do 
not enable to perform clustering measurement separately 
on these populations. 

We computed the clustering of galaxies within the 
green valley as defined in sect. 3] and using dust correc- 
tion, but separately for galaxies with or without AGNs. 
The AGNs represent 56% of the green valley. The pro- 
jected correlation functions are presented in fig. 1131 The 
main differences in terms of clus tering between th ese two 
population, as noted already bv lLi et all (|2006bD . are at 
small scales. At scales 0.1 < r v < 0.5/z." 1 Mpc, the non- 
AGNs galaxies exhibit marginally more pairs than the 
AGNs. At large scales, the clustering of AGNs and non- 
AGNs populations is very similar; the correlation length 
for the AGNs is 4.7 ± 0.3/i _1 Mpc, and the correlation 
length for the non-AGNs is 5.1 ± 0.4/i _1 Mpc. This sug- 
gests that the intermediate clustering of the green valley 
galaxies is genuine and not linked to the predominance of 
AGNs in this region of the UV optical color-magnitude 
diagram. 

5.2. Dark Matter Halo mass 



We used the analytical models of iMo fc Whitel (|2002l ) 
to estimate the minimum mass of the dark matter 
halos that have similar clustering, by matching the 
bias parameter fr 8 , computed at 8 /i _1 Mpc, following 
e.g. iMagliocchetti et al.l (|2000f ). This comparison shows 
that the mass of the host dark matter halos increases 
smoothly with (NUV — r) a b s - Blue galaxies have clus- 
tering similar to halos with minimum masses around 
lO n /i _1 M0, green valley galaxies 10 12 - 5 ft. _1 M Q , and red 
galaxies 10 13 /i -1 Mq. The minimum halo mass for blue 
galaxies is in agreement with estim ates derived from 
UV-selected galaxies at low redshift (jHeinis et al.l 120071 
Millia rd" et alj l2007f) . and is additional evidence that 
active star formation preferentially occurs in low den- 
sity environment s at low redshifts (| Gomez "et~all[200l 
ILewis et al1l2002D . 

The minimum dark halo mass for transition galaxies 
is similar to those of galaxy groups. This suggests that 
galaxies within the green valley mainly reside in groups. 
At z ~ 0.85 gr een valley galax ies have a higher bias, 
b = 1.65 ± 0.26 (|Coil et al J 12008): at these redshifts this 
corresponds also to halos with minimum masses of the 
order of 10 12 ' 5 /i _1 Mq. The combination of these results 
show that green valley galaxies reside in groups of similar 
masses at z ~ 0.85 and z = 0. We note also that within 
the green valley the correlation function of star forming 
galaxies is marginally larger at small scales than the AGN 
correlation function (see fig. [13]) . This could provide 
some hints about the nature of those objects in the con- 
text of the Halo Occupation Distribution and the satellite 
versus central galaxy population. As AGNs are expe cted 
to be preferentially central galaxies (|Li et al.ll2006b[ ) this 
is a hint that transition galaxies could be mainly satel- 
lites rather than central galaxies within groups and may 
preferentially lie in the outskirts of their host DMHs as a 
recent ly accreted population, as proposed by iCoil et al.l 
(2008). We can speculate on the nature of the mecha- 
nisms besides AGN feedback involved in the quenching of 
star formation in transition galaxies, if they are related 
to their environment. In this case, given the low relative 
velocities in galaxy groups, local processes such as tidal 
mechanisms m ay induce significant suppression of star 
formation (e.g. Bosclli & Gavazzi 2006; Bowc r"fc Baloghl 
l200llG6mez et all 12001 ILewis et al.ll2002l) . 

Finally, as expected, the reddest galaxies have a clus- 
tering strength similar to large groups or clusters. 

These results show that there is a correlation be- 
tween star formation history and dark matter halo mass, 
as the galaxies having formed stars more recently in- 
habit less massive halos. These fin dings are in agree - 
ment with the "downsizing" scenario (|Cowie et alJl 996). 
which states that more massive ga laxies formed their 
stars at higher red shifts (see e.g. iBundv et al.l 120061 : 
iDe Lucia etaTll2006l ). 

5.3. Halo occupation limits 

Our approach in the previous section to estimate 
the dark ma t ter h alo mass from the clustering with 
IMo fc White! (|2002l ) models has the caveat that it as- 
sumes only one galaxy per halo. This assumption is 
quite unlikely accordi ng to recent numerical and obser- 
vatio nal studies (e.g.. iKravtsov et "all 120041 : lOuchi et al.l 
2005). In the frame work of the Halo Occu pation Distri- 
bution (HOD) (e.g.. lCoorav fc Sheth|[2002T ). the number 
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Fig. 14. — Bias as a function of galaxy number density for 
(NUV — r) a bs cut samples (same color legend as fig. [HJ. We 
show here only results without the dust extinction correction. 
The solid li ne shows the relatio n expected for dark matter 
halos, using iMo fc White! (|2002t l models. 

of galaxies i n a given halo depen ds on its mass. We can 
however use IMo fc White! (|2002f ) models to get a first 
estimate of the dependence of the halo occupation with 
star formation history. In fig. [14] we compare the bias 
and the number density of the samples with those of the 
dark matter halos. We computed the number densities 
of the v arious samples using the luminosity functions de- 
rived bv lWvder et all {2007) in (NUV - r) abs bins from 
a similar MIS sample. These luminosity functions are 
not available for dust extinction corrected samples, so 
we only show on fig. [14] our MIS results obtained with- 
out any dust correction. As the clustering results are in 
very good agreement with or without applying the dust 
correction, and as we want to get only some qualitative 
limits on the halo occupation numbers, we expect our 
conclusions to be robust. 

The relation expe cted for dark matter halos from the 
IMo fc White! (|2002D models is shown by the solid line. 
On this plot, the dark matter halo mass increases on this 
line from right to left. As this model assumes that each 
halo hosts one galaxy, measurements on the solid line 
indicate an average halo occupation number of 1, mea- 
surements above an average halo occupation greater than 
1, and measurements below an average occupation lower 
than 1. Our measurements recover the trend expected 
for dark matter halos, as galaxy populations with lower 
number densities are more clustered. However, galaxies 
do not lie on the dark matter halos bias-density rela- 
tion, which is due to the fact that the halo occupation 
varies among populations. Our results show that the av- 
erage occupation number increases with (NUV — r) a b s 
color. Bluest galaxies in the local Universe have a very 
low occupation number in halos that have similar clus- 
tering strength. Our second blue sample contains mainly 
galaxies with an average occupation number close to 1. 
Transition and red galaxies have an occupation number 
around 5. Note that the occupation number for the red 
galaxies is potentially an under estimation given that our 
selection biases against low mass and low specific star 



Fig. 15. — Specific star formation rate as a function of stellar 
mass. The dot-dashed contours show the distribution of the 
green valley defined with 2.75 < (NUV — r) a b s ,corr < 3.5. 
The dashed lines represent the cuts we consider for clustering 
measurements. The error bars show the dispersion caused by 
the dust extinction correction, 
formation rates galaxies (see sect. I4.2[) 

These r esults are in qualit ati ve agreemen t s wit h 
numeric (jBerlind et all 120031 : IZheng et all I2005T ). 
as well as observational studies (iHogg et al.l 120031 : 
iMagliocchetti fc Porcianill2003l : IZehavi et al.ll200'oT) . 

5.4. Specific star formation rate vs. stellar mass 

The results we present above add further evidence 
that the clustering of galaxies strongly depends on star 
formation history. Howe ver galaxy clus t ering also de- 
pend s on stellar mass (e.g. Hi et al.|[2006al : iMeneux et all 
2008); we investigate which of these two parameters 
shows the strongest correlation with environment. To do 
so, we consider an alternative to the UV-optical color- 
magnitude diagram, by plotting the specific star forma- 
tion rate SFR/M* as a function of the stellar mass M* 
in figure [TS] 

We computed the star formation rates from the NUV 
lumin osities correcte d for dust extinction (see P2.3"]) using 
the lKennicuttl (| 19981 ) relation scaled to a lKroupal (|200ll ) 
initial mass function: 



SFR (A/qvi'- 1 ) = 10- 28 02 L„(erg s" 1 Hz" 1 ). (8) 

We est imate the st e llar m asses using the equation de- 
rived by lYang et"afl (|2007l ) from the relation betwee n 
stellar mass-to-light ratio and color of iBell et~ al. ( 2003): 



log 



M, 



-0.306 + 1.097[(g o6 



-0.1 



Tabs ) 



QA[r abs ,oorr - 4.64] (9) 

which implicitly assumes a lKroupal (|2001l ) initial mass 
function. 

The cuts we adopted are shown as dashed lines in fig. 
IT5l Within the stellar mass range 10 9 9 < M* < 10 12 Mq, 
we choose an additional mass cut at M* = 3 x 1O 1O M . 
Several studies indeed showed that in the local Universe 
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TABLE 4 

MIS/SDSS sample M» — SSFR clustering results* 





Sample + 






Ngal 


<*> 


<M,)[M ] 


r [h -1 Mpc] 


7 


b 8 






10 8 - 


< M, < 10 9 9 


2103 


0.04 


10 9 56 


2.46 ±0.21 


1.6 ± 0.1 


0.56 ± 0.04 


-10. 


< \og(SFR/M t ) < -9. 


10 99 


< M, < 10 10 ' 48 


3347 


0.07 


10 10.24 


3.58 ±0.17 


1.65 ± 0.07 


0.76 ± 0.03 






^glO.48 


< M, < W 12 - 


3375 


0.12 


^ gl0.82 


3.48 ±0.16 


1.65 ± 0.06 


0.76 ± 0.03 


-11. 


< log(SFR/M*) < -10. 


10 99 
1Q10.48 


< M* < 10 1048 

< M, < 10 12 ' 
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Fig. 16. — Projected correlation functions of the MIS/SDSS 
samples cut by stellar mass and specific star formation 
rate. The top panel shows measurements obtained for the 
-12.6 < log(SSFR) < -11 bin, the middle panel for the 

— 11 < log(SSFR) < —10 bin, and the bottom panel for the 

— 10 < log(SSFR) < —9 bin. On each panel, filled squares 
show results from 9.9 < log(M») < 10.48, filled circles from 
10.48 < log(Af„) < 12, and filled triangles 8 < log(M») < 9.9 
(only for -10 < log(SSFR) < -9). 

this value divides galax ies between active or pa ssive in 
terms of star formation ( Kauffm ann et al. 2003a) . More- 
over, according to numerical studies ( Keres et al.ll2005i 
2008), the way galaxies accrete gas depends on stellar 
mass. Galaxies below 3 x 10 10 M Q are dominated by a 
"cold mode" of gas accretion, while galaxies above this 
mass are dominated by a "hot mode" of gas accretion. 
The geometry of these two modes are also different, the 
cold mode being mostly filamentary, while the hot mode 
quasi-spherical, hence implying a link with environment. 
We then further divide our subsample in mass by this 
value M* = 3 x 1O 1O M0, and also use two cuts in specific 
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Fig. 17. — Correlation length as a function of specific star 
formation rate and stellar mass. The filled triangle shows the 
result from 8. < log(M») < 9.9, filled squares show results 
from 9.9 < log(M„) < 10.48, and filled circles from 10.48 < 
log(M„) < 12.. 

star formation rate at \og(SFR/M 1 ,[yr^ 1 ]) = 10~ n and 
logiSFR/M^yr- 1 }) = 10~ 10 . 

The correlation functions are shown fig. [TH] Each 
panel of fig. [TB] shows a given cut in specific star for- 
mation rate, and the relevant cuts in stellar mass. While 
there are some differences in the small-scale clustering 
(r p < 0.5/i _1 Mpc), at larger scales the correlation func- 
tions for the galaxies with 10 9 9 < M* < 10 12 M Q are re- 
markably similar for a given specific star formation rate. 
The main difference we observe in clustering strength 
is for the less massive sample (10 8, < M„ < 10 9 ' 9 M Q ) 
within the highest specific star formation rate bin. There 
is also a difference for the lowest specific star forma- 
tion rate bin, but in this bin the low mass sample with 
10 9 - 9 < M* < 10 10 - 48 M Q has very low statistics so its 
correlation function must be interpreted with caution. 

We show in fig. [IT] the correlation lengths derived 
from these measurements as a function of specific star 
formation rate and stellar mass. The global trend is 
that rn decreases with specific star formation rate. For 
-11 < ]og(SFR/M») < -9 and 9.9 < log{M*) < 12, 
the correlation length does not depend on stellar mass 
in a given specific star formation rate range. This does 
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not seem to be true in our extreme stellar mass bins. 
However the statistics of our -12.6 < log(SFR/M*) < 
— 11, 9.9 < log(M») < 10.48 sample are small; and 
because of our selection, we can not compare to lower 
specific star formation rates at 8. < log(M») < 9.9. Nev- 
ertheless, these results suggest that the specific star for- 
mation rate shows a stronger correlation with environ- 
ment (as traced by the large-scale clustering) than the 
stellar mass. Interestingly, this is true below and above 
the transition mass of 3 x 10 10 Mq. Within the ranges 
-11 < \og(SFR/M*) < -9 and 9.9 < log(M*) < 12, 
our results show that there is a smooth transition of the 
large-scale clustering strength, and hence in dark matter 
host halo masses, as a function of star formation history, 
independently of stellar mass. This result is consistent 
with our measurements based on color selected samples 
(because of the correlation between specific star forma- 
tion rate and (NUV — r) a b s ), but provides a direct in- 
terpretation of the dependence of clustering on galaxies 
physical parameters (stellar mass and specific star for- 
mation rate). 

In the local Universe, high specific star formation rate 
galaxies are anti-biased (6§ < 1), and have clustering sim- 
ilar to dark matter halos of minimum masses ~ 10 11 M Q . 
Galaxies with intermediate specific star formation rates 
and within the green valley are unbiased (bg ~ 1), and 
arc hosted by dark matter halos of minimum masses 
~ 10 12 Mq. Lower specific star formation rate galaxies 
are biased (6s > 1); their clustering strength is similar to 
dark matter halos of minimum masses ~ 1O 13 M . 

These resu lts are in agreements with those of 

iBlanton et al.l (|2005f ) who noticed that color (g — r in 
their case, which is also a tracer of star formation his- 
tory) i s the galaxy prop e rty m ost predictive of environ- 
ment. iKauffmann et all (|2004f ) showed that the galaxy 
property most sensitive to environment is star forma- 
tion history (and in particular the specific star forma- 
tion rate), and although the amplitude of the depen- 
dence varies with mass, it is observed at all masses (for 
10 9 - 5 < M* < 10 n M Q ). 

6. CONCLUSIONS 

We used the unique combination of GALEX UV pho- 
tometry and SDSS spectroscopy to perform spatial clus- 
tering studies from UV and UV/optical-selected samples. 
We take advantage of the remarkable correlation between 
the (NUV — r) a bs color with star formation history to in- 
vestigate directly how it relates to large-scale clustering. 

The results fro m the UV-selected sa mples, in addi- 
tion to those from lHeinis et all (|2007f ) and lMilliard et all 
(2007), confirm the low clustering of UV-selected galaxies 
in the local Universe. This is further evidence that ac- 
tive star formation takes preferential ly place in low den- 
sity environment s at low redshifts (j Gomez etHl [20031 
iLewis et~al"1l2002j ). Within the blue sequence, redder UV- 
selected galaxies are more clustered, but the clustering 
strength of UV-selected galaxies is however lower than 
for optically-selected galaxies at a given (NUV — r) a b s 
color. 

The results from our MIS/SDSS sample recover the ex- 



pected trends from previous studies that redder galaxies 
are more clustered. Thanks to the greater sensitivity of 
the (NUV — r) a bs color to star formation history than 
optical colors, our results span a wide range in clustering 
strength. Our results show that galaxies that experi- 
enced star formation more recently inhabit less massive 
halo s, which is consistent with the "downsizing" scenario 
(e.g.lBundv et al.ll2006tlCowie et al.lll996l : lNeistein et all 



2006). Also, their average halo occupation is lower, in 
such a way that the bluest galaxies are hosted by only 
a small fraction of the dark matter halos that have the 
same clustering. These findings show that, in contrast 
to what is observed from LGBs at high redshi ft (see e.g. 
iGiavalisco fc Dickinson! 120011 : iLee et all 12006). star for- 
mation is not primarily driven by gravitation in the lo- 
cal Universe. The clustering of the green valley galaxies 
suggests that they are mainly satellites wit hin galaxy 
group s. Comparison with previous work of ICoil et all 
(2008) indicate that these transition galaxies reside in 
similar environments at z ~ 1 and z — 0. 

Our results show that there are smooth trends between 
star formation history and environment. It is interesting 
to notice that based on the results from optical selections, 
one might not expect to observe such smooth transitions. 
Previous studies showed that there is a bimodality in a 
nu mber of galaxy properties: color (u — r for instance; 
seelBaldrv et al.ll2004t IStrateva et al.ll200lb. stellar mass 
(IKauffmann et aLl2003bh. or clustering (jBudavari et al.l 
120031 ). IWvder et all ( |2007t ) already noticed however that 
the (NUV — r) a bs color distribution in the local Universe 
is not as bimodal as the u~r color (see also our fig.H]), due 
to the greater sensitivity of the UV spectral range to star 
formation activity, and reddening by interstellar dust. 
Our results suggest that in the local Universe, rather 
than a bimodality, there is a smooth transition in global 
environments (as probed by the clustering) for galaxies 
with different star formation histories (as probed by the 
(NUV — r) a bs color). These results hold whether dust 
extinction is applied or not, and are insensitive to the 
predominance of AGNs among transition galaxies. 

Finally, our results show that the large-scale clustering 
trends we observe are mainly driven by star formation 
history, and not stellar mass. This is indication that at a 
fixed star formation history, there is no strong variation 
of the host dark matter halo mass. Our results suggest 
that the environmental dependance is first driven by the 
star formation history and that the stellar mass only acts 
as a secondary parameter. 
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